Droplets of mixtures of 7He-*He at finite vorticity 
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Abstract 


Stimulated by recent interests in mixtures of quantum fluids, we study the properties of He-+He mixture droplets 
using a phenomenological density functional. The structure of liquid drop with and without vortices is investigated 


using Monte Carlo simulations. 
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There is a burgeoning interest in droplets of quan- 
tum fluids from different directions recently. These in- 
clude droplets of Bose-Einstein condensed alkali atoms 
confined by magnetic and optical means[1], small 
droplets (1000 atoms) of He in a molecular beam[2], 
and large droplets (mm in diameter) levitated by a 
magnetic field[3]. While one component drops have 
received much interest, multicomponent systems are 
beginning to be looked at. In addition to their own in- 
trinsic interest, droplets of mixtures provide a unique 
opportunity to test and enhance our understanding of 
finite size effects and the influence of boundaries on 
phase segregation effects. 

Current Bose condensed alkali mixtures are mostly 
mixtures of Bose particles. 7He-*He are mixtures of 
Fermions and Bosons. In this work we present Monte 
Carlo simulations on the structure of 7He-*He mixture 
droplets with and without vortices as a function of their 
relative concentration. This follows our earlier work 
that focuses on the simpler aspects of this system|[4]. At 
zero vorticity, because of the quantum pressure most 
of the *He stays on the outside with a small amount of 
3He remaining inside, as is previously reported. 

The total free energy for a liquid-helium mixture of 
3He and *He can be expressed in a density functional 
form.[6] We approximate the volume integral of the 
energy functional by a discrete sum. We sample a lat- 
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Fig. 1. Cross sections of density distributions of 3He in the xy 
plane for z=0 with particle numbers N3 = 500 and N4 = 3500 
(dashed line), N3 = 3000 and N4 = 1000 (solid line) 


tice inside a sphere of radius 2R consisting of 40 sites 
along the diameter, making a total of 33398 sites. We 
choose the mesh size as 3 A, which is smaller than the 
surface and interface thickness. We use the Metropo- 
lis algorithm, start our simulation with high tempera- 
ture, then decrease it to a very low temperature step 
by step. This is an annealing process and the final con- 
figuration should be the lowest energy configuration. 
Because of the strong interactions between the helium 
atoms, the particle number is a very sensitive function 
of the chemical potential. We thus carry out the simu- 
lations with fixed particle numbers. 
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Fig. 2. Image of density distributions of ?He with a vortex 
in “He. The vorticity is m = 1 and the particle numbers are: 
N3 = 3500 and N4 = 500. 


Cross sections of the snap shots of the density pro- 
file for two different concentrations are shown in Fig. 
1. The maximum bulk solubility of *He in the “He-rich 
phase is ~ 6.6% at zero temperature. When the 7He 
fraction is small, there is a small amount He inside 
the droplet, consistent with the bulk phase diagram. 
However, the concentration of *He is less than that in 
the bulk. We think this comes about because the *He 
inside must have the same chemical potential as that 
outside. Because of the free boundary condition the 
chemical potential of “He outside is lower than that in 
the bulk. When the *He number is very small, the sol- 
ubility of He into “He is much smaller than the bulk 
value. The effective attractive potential for 7He atoms 
on the surface of “He provide a type of surface states, 
which are called Andreev states. “He atoms first accu- 
mulate on the free surface occupying Andreev states 
before start dissolving into the bulk. The recent change 
of the He3-4 mutual solubility in aerogels[5] is another 
example of the effect of boundary conditions on the 
phase segregation phenomena. 

Taking the vortex along the z axis, the free energy 
functional of the vortex state is obtained with the ad- 
dition of a centrifugal term A? pak? /(2mar’) associated 
with the superfluid flow where « is the quantum circu- 
lation number and r is the distance to the z axis. At 
finite vorticity, “He moves outside because of the cen- 
trifugal force. Images of the snap shots of the density 
distributions at a low “He concentration (N3 = 3500 
and N4 = 500) are shown in Fig. 2. The “He forms 
lumps and does seems to “wet” the He completely. We 
emphasize that the pictures are snap shots and not 
averages of the denisty distribution. It is the density 
correlation between He and “He that is being shown. 
Thus at low “He concentrations, the density correla- 
tion between *He and “He becomes asymmetric. Since 


“He is rotating, we expect the average density to be 
cylindrically symmetric. We think this phenomenon is 
related to the case of bulk 7He-*He mixture, where the 
votices lose their stability with increasing He concen- 
trations[?,?]. The more interesting question is how big 
will the damping be against the superfluid flow and 
how metastable is this configuration. We hope to in- 
vestigate this in future publications. 
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